CGC predictions for p+A collisions at the LHC and signature of QCD saturation 
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We present various predictions for the upcoming p+Pb collisions at \/S = 5 TeV within the color 
glass condensate (CGC) formalism, including single inclusive charged hadron production, single 
inclusive prompt photon production, direct photon production, charged hadron multiplicity distri- 
bution and photon-hadron azimuthal correlations. Using the running-coupling Balitsky-Kovchegov 
evolution equation for calculating various observables, we show that the main source of uncertainties 
is due to less constrained initial nuclear saturation scale. This gives rise to rather large theoretical 
. . . uncertainties for nuclear modification factor RpA at the LHC. Nevertheless, we propose a simple 

(N- scheme in which one can still test the main dynamics of the CGC/saturation in p-l-A collisions at 

^ ■ the LHC. 

o: 

(N. 

^ ■ I. INTRODUCTION 

The upcoming proton-lead (p-|-Pb) collisions at the LHC will provide crucial benchmarks for understanding the 
characteristic of Quark Gluon Plasma produced in heavy ion collisions at the LHC and RHIC. Moreover, p-l-A collisions 
\ run has its own merits, namely it can be vital testing grounds for novel nontrivial QCD dynamics which otherwise it 
cannot be unambiguously explored. It is generally believed that a system of partons (gluons) at high energy (or small 
Bjorken-x) forms a new state of matter where the gluon distribution saturates and non-linear coherence phenomena 
dominate Such a system is endowed with a new dynamical momentum scale, the so-called saturation scale 
which controls the main characteristic of the particle production. The color glass condensate (CGC) approach was 
proposed as a framework to study small-x and saturation physics [5]. The CGC formalism is an effective perturbative 
^ ■ (weak-coupling) QCD theory in which one systematically re-sums quantum corrections which are enhanced by large 
logarithms of 1/x and also incorporates high gluon density effects which are important at small x and for large nuclei, 
for a review see e.g. ReL [1]. 

(Nl I The CGC formalism has been successfully applied to many processes in high energy collisions. Examples are 
structure functions (inclusive and diffractive) in Deeply Inelastic Scattering of electrons on protons or nuclei, and 
particle production in proton-proton, proton-nucleus and nucleus-nucleus collisions at RHIC, see Ref. Q and references 
. therein. The observed suppression of the single inclusive hadron production in deuteron-nucleus collisions (d-|-A) at 

■ forward rapidities at RHIC [3] is perhaps among the most spectacular evidence of the QCD gluon saturation. There 
, are, however, alternative phenomenological approaches which describe the same data. Therefore, the upcoming p-l-A 

run at the LHC which covers a much wider kinematic region, can provide a unique opportunity to understand the 
I underlying dynamics of particle production at small x. 

The key ingredient of the particle production at leading logarithmic accuracy which captures the main saturation 
dynamics is the universal color dipole cross-section, the imaginary part of the quark-antiquark scatterings amplitude. 
It is universal, in a sense that for different processes, the cross-section can be written in terms of the same object, 

■ see e.g. Refs. @, @|- It has been recently proven that in the large- A'c limit, all multi-particle production processes 
' in the collision of a dilute system off a dense (e.g, p-l-A collisions) can, up to all orders in the strong-coupling, be 

described in terms of only dipoles and quadrupoles (quadrupoles can be also approximated to dipoles) [6|. The 
color-dipole amplitude satisfies the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution 
equations H . In the large Nc limit, the coupled JIMWLK equations are simplified to the Balitsky-Kovchegov (BK) 
equation [a Q, a closed- form equation for the rapidity evolution of the dipole amplitude. Unfortunately, at the 
moment, we do not know how to deal with impact-parameter dependence of the BK (or JIMWLK) evolution equation 
[Tol [Tl|. and impact-parameter dependence of the dipole amplitude is generally ignored. Another current drawback 
is that numerical solutions of the full next-to-leading logarithmic expressions [l3| are not yet available, only running 
coupling corrections to the lea ding log kernel have been considered in phenomenological applications, the so-called 
running-coupling BK equation \V^. Moreover, the BK evolution only provides the rapidity/energy evolution of the 
dipole, the initial profile of the dipole amplitude and its parameters still need to be modeled and constrained by 
experimental data. Unfortunately, the current world-wide small-x data are very limited, and cannot uniquely fix the 
initial parameters of the dipole amplitude . This problem is more severe for determining dipole scattering amplitude 
on nuclear target. This leads to rather large unavoidable theoretical uncertainties for various CGC predictions in p+A 
collisions at the LHC. Here, we propose a simple scheme in which one can overcome this problem and still test the 
model at the LHC. 

In this letter, we provide a compilation of various predictions for p+Pb collisions at ^/S = 5 TeV at different 
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rapidities based on the CGC approach. This includes the total charged hadron multiplicity distribution at different 
centralities, nuclear modification factor for single inclusive charged hadron production, nuclear modification factor for 
single inclusive prompt photon production (and direct photon production), and azimuthal photon- hadron correlations. 
These observables have been already studied in the CGC framework but at different energies at the LHC, namely 
= 4.4 and 8.8 TeV [IMS- 

Here we also extend the previous studies by quantifying theoretical uncertainties 
coming from our freedom to choose different values for the factorization scale Q, the strong-coupling as, and initial 
saturation scale of proton and nucleus. 



II. THEORETICAL FRAMEWORK AND MAIN FORMALISM 

A. Single inclusive hadron production in p+A collisions; kr factorization, and hybrid formalism 

In the CGC approach, the gluon jet production in p-l-A collisions can be described by /cT-factorization [l9l|, 

da 2q!, 1 



dytPpT Cp Pt 



d^kT(t>p {xi;kT^ (j)'^ (x2;pt - ^t) , (1) 



where Cp = (N^ — ^)/'2,Nc with Nc being the number of colors, xi^2 — {PT/VS)e^y, pT and y are the transverse- 
momentum and rapidity of the produced gluon jet with \/S being the nucleon-nucleon center of mass energy. 
(ff^{xi;kT) denotes the unintegrated gluon density and is the probability to find a gluon that carries Xi fraction 
of energy with fcr transverse momentum in the projectile (or target) A. The unintegrated gluon density is related to 
the color dipole forward scattering amplitude. 



05 (x,- kr) = ^ / d^br d^rr^^^^-^^^lUA [x,; tt'. br) , (2) 

with notation 



J\fA (xi^rT^br) = 27V> {xi;rT;bT) ~ Afp {xi\rT\bT) , (3) 

where tt denotes the dipole transverse size and br is the impact parameter of the scattering. Throughout this paper, 
the subscript T stands for the transverse component. For the value of strong-coupling as we employ the running 
coupling prescription used in Refs. j20l - [23j . Namely, in Eq. (jlj we replace as by agipx), and in Eq. ([2|) we replace as 
by as{Qs{xi)) with Qs{xi) being the saturation scale in the projectile or target j20i] . The most important ingredient 
of the fcr factorization which encodes the saturation dynamics is the fundamental (or adjoint) dipole amplitude, the 
imaginary part of the forward quark anti-quark scattering amplitude on a proton or nucleus target Np^A {xi]rT]b). 
The dipole cross-section can be computed via the perturbative nonlinear small-x Balitsky-Kovchegov (BK) quantum 
evolution equation see below. Note that the /ct factorization has been proven at the leading log approximation 
for scatterings of a dilute system on a dense one (such as proton-nucleus collisions) and includes BFKL type gluon 
emissions with gluon fusion effects between the projectile and target, and also gluon radiations from the produced 
gluons flQ*!. This formulation has been independently verified by many authors, see e.g. Ref. [23 |. 
The rapidity distribution of the inclusive mini-jet production can be calculated from Eq. ([T|) , 

dN\„^ K f dn 

jet _ _ / d^p^h[^^p^^^^^^^^^[Eq. ©], (4) 



dr] Us J ' ' dy d^pT 

where "q is the pseudorapidity and h\r), pj, rnne.t] is the Jacobian which takes account of the difference between rapidity 
y and the measured pseudorapidity r] [20|. In the above, we introduced a factor to effectively incorporate the 
missing higher-order corrections, as is the average interaction area which depends on the kinematics and is related 
to the average impact parameter of the inclusive production of the mini-jet ^20,] . By introducing mini-jet mass uijet 
which mimics the pre-hadronization effect, one can regularize the infrared divergences of the kx factorization cross- 
section [l^ [20I - I23I I. Finally, in order to relate the mini-jet yield to the produced hadron, we employ the so-called 
Local Parton-Hadron Duality principle ^25] assuming that the form of the rapidity distribution will not be distorted 
by the jet decay and only a numerical factor will differ the mini-jet spectrum from the hadron one. Note that the 
main contribution of the integrand in Eq. Q comes from low transverse momentum (typically less than 1-^3 GeV) , 
therefore using fragmentation function is less justifiable. 

In the fcj^-factorized approach, partons in both the projectile and target are assumed to be at very small x in order 
the CGC formalism and small-x resummation to be applicable to both the projectile and target. This approach is 
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valid as long as one stays away from the projectile fragmentation re gion . However, for the projectile fragmentation 
region at forward collisions, an alternative approach was developed in [26ll27l[. the so-called hybrid approach, which is 
better suited where one treats the projectile wave-function perturbatively within the standard coUinear factorization 
approach using the standard DGLAP picture while treating the target by CGC methods. The cross section for single 
inclusive hadron production at leading twist approximation, in asymmetric collisions such as p-|-A ones, in the CGC 
formalism is given by [IB, US] , 



(P-prdri ^ (27r)2 
a 



1 dz 

XF 

1 dz Z^ 



lfg{xi,Q )NA{x2,—)Dh/g{z,Q) + T.qXlfq{xi,Q )NF{x2,—)Dh/q{z,Q) 

/ d^kTkpNF{kT,X2) / —Y:i^j=q,q,gW^/j{^)P,_/j{C)Xifj{-^,Q)Dh/i{z,Q) 
Jk^<Q^ Jxi t. t. 



,(5) 



where a™ behind the inelastic term is the QCD strong-coupling and denoted by subscript "m" in order to be 
different with the strong-coupling as-running in the rcBK equation. Note that in the hybrid formulation given above, 
in principle the strong coupling in the projectile which is in dilute regime can be different from the one in the rcBK 
description of dense target. The scale of a™ cannot be determined at its current approximation, and for that a full 
calculation upto NNLO is required. We will later consider the implication of various values for a™. fj{x,Q^) is 
the parton distribution function (PDF) of the incoming proton which depends on the light-cone momentum fractions 
X and the hard scale Q. The function Dfi/i{z,Q) is the hadron fragmentation function (FF) of i'th parton to the 
final hadron h with a momentum fraction z. The indices q and g denote quarks and giuon, with a summation over 
different fiavors being implicit. The inelastic weight function Wi/j and the DGLAP splitting functions Pi/j are given 
in Ref. (27j . The longitudinal momentum fractions xi and X2 are defined as follows. 



Xp 



Pt 



X\ 



Xp 

z 



X2 



-2r, 



(6) 



Note that in the hybrid formalism Eq. ([5]), in contrast to the kr factorization Eq. ([TJ, hadronization can be treated 
similar to the standard pQCD due to resummation of collinear singularities via DGLAP evolution for the incoming 
and outgoing parton. The expersion given in Eq. ([5]) can be understood in a simple intuitive picture. The first two 
terms correspond to elastic contribution, namely an incoming parton scatters elastically with the CGC target. This 
incoming parton with initial zero transverse momentum picks up transverse momentum of order saturation scale after 
multiple scatterings on the dense target. The last term in Eq. ([5]) gives the inelastic contribution to the inclusive 
hadron production. Namely, the projectile parton can also interact with target inelastically with small transfer 
momentum exchanges. This decoheres the pre-existing high-pr parton from the hadron wave function and releases 
it as an on-shell particle. These high-pT partons in the projectile wave function arise due to DGLAP splitting of 
partons. 

In above, the amplitude Np (Na) is the two-dimensional Fourier transformed of the imaginary part of the forward 
dipole-target scattering amplitude AfA(F) in fundamental (F) or adjoint (A) representation. 



NAiPj{x,kT) = / d^fe-'^^-''{l-MAiP){r,Y = Hxo/x))) , 



(7) 



where r = is the dipole transverse size. The dipole forward scattering amplitude incorporates small-x dynamics 
and can be calculated via the running-coupling BK (rcBK) evolution equation. The rcBK equation has the following 
simple form: 



= / d^fi K'-'''^{r,ruf2)[J^A{F){rux)+NA{P){r2,x)-JVA(P)ir,x)~~UA(F){ri,x)J\fA^ (8) 



dAfA(P){r,x) 
dln{xo/x) 

with r2 = r — ri- The only external input for the rcBK non-linear equation is the initial condition for the evolution 
which is taken to have the following form motivated by McLerran-Venugopalan (MV) model [2,] , 



7V(r,r = 0) = 1 - cxp 



2 r)2 \T 



In ( + e 
Ar 



(9) 



where the onset of small-x evolution is assumed to be at xq = 0.01, and the infrared scale is taken A — 0.241 GeV 
14| . The only free parameters in the above are 7 and the initial saturation scale Qqs, with a notation s = p and, A 
for a proton and nuclear target, respectively. We will later consider uncertainties coming from our freedom to choose 
among different parameter sets of the rcBK description of the proton and nucleus. 
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B. Direct photon production and photon-hadron correlation in p+A collisions 

The cross section for semi-inclusive prom pt p hoton-quark production in p-|-A coUisions at the leading twist approx- 
imation in the CGC formalism is given by (28l [29|. 



(Pbr (PpV <PIt dvj dr^h V2(47r4) [p^f^ [p- 1^ - l-pV? 

2rp- It ■ pV + {k- -p-)lT + l'{k- - H (p^f 



Npilh+PT^lxg), (10) 



where p'^ , I, and k are 4- momenta of the produced prompt photon, outgoing quark and projectile quark, respectively. 
Again, a ii'-factor was introduced to absorb higher-order corrections. The light-cone fraction Xq is the ratio of the 
incoming quark to proton energies, namely Xq = k~ j \J S/2. The pseudorapidities of outgoing prompt photon 77-y 
and quark 77^ are defined via = ^^''^ ^^'^ ~ '^^^ angle between the final-state quark and prompt 

photon is denoted by and defined via cos(A0) = 7^-^- We only consider here light hadron production, therefore 

at high transverse momentum, the rapidity and pseudorapidity are the same. In Eq. pH)) . Np{pT,Xg) is again 
Fourier transformed of the dipole amplitude which satisfies the rcBK equation The semi- inclusive photon-hadron 
production in proton-nucleus collisions can be obtained from partonic cross-section Eq. (|10p by convolution of quark 
and antiquark distribution functions of a proton and the quark-hadron fragmentation function, 

= Z ~ = / / dXqfq{Xq,Q^) Dh/q{zf,Q^), (11) 

d^hT d^pr'' d^PT drj^drjh J zf^^ ^/ J dPbx d'^pr'' cPIt drj^ dr]h 

where p^ is the transverse momentum of the produced hadron. A summation over the quark and antiquark flavors in 
the above expression should be understood. The light-cone momentum fraction Xq,Xq,Xg in Eqs. (|10llip are related 
to the transverse momenta and rapidities of the produced hadron and prompt photon via fl7| , 

^ f ^ V ^ 

Xq — Xq — -p=. I prp e ^ H e ' 

^ ( ^1 ^-v-, ,Pt_ ^-Vh 



Zf 



V Zf 

P%Ht, with = ^ I 1 . (12) 



V5 



7s 



The single inclusive prompt photon cross section in the CGC framework can be obtained from Eq. (jlOp by integrating 
over the momenta of the final state quark. The cross-section of single inclusive prompt photon production can be 
divided into two contributions of fragmentation and direct photon jl7| : 



^^qA— ^7(p^)X ^^Fragmentation ^^Dircct 



d'^brd^p'r'^ drij d^brd'^PT^ drj-f d^bTd^pr^ d-q-^ 
K 



(13) 



1^ / ^Tv,^/ -vys ^n^era 



D,,q{z,Q^)NF{Xq,pl/z) + ^i-^Z^[l + {l-zf]y^^ / ;2 ( ^ ) 



where D^/q(z, Q^) is the leading order quark-photon fragmentation function 30]. Similar to the hybrid formalism for 
the inclusive hadron production Eq. ([5]), Q is a hard-scale. One should note that above expersion was obtained using 
a hard cutoff to subtract the collinear singularity [Tt'I . This may result in a mismatch between the finite corrections to 
our results and those that are included in parameterizations of photon fragmentation function. However, this mismatch 
is a higher order effect in the coupling constant and its proper treatment requires a full NLO calculation which is 
beyond the scope of this letter. In order to somehow quantify possible errors associated with the approximation made 
in Eq. (jl3p . we will consider different photon fragmentation functions, and also different values for the hard-scale Q. 
In order to relate the partonic cross-section given by Eq. (jl3p to prompt photon production in p+A collisions, we 
convolute Eq. (|13p with quark and antiquark distribution functions of the projectile proton, 

Z = / dXqfq{Xq,Q^) , (14) 

d^bTd^Pr^ drj^ Jx'^^^ d^bxd'^Px'^drij 
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where a summation over different quarks (antiquarks) flavors is implicit. The light-cone fraction variables Xg,Xg,z in 
Eq. (jl3ll4p are related to the transverse momentum of the produced prompt photon and its rapidity {T^ . 
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Pt 

7^' 



(15) 



Similar to the inclusive hadron production, the main ingredient of photon-hadron and prompt photon production 
cross-section is the universal dipole amplitude. One should note that the light-cone fraction variables which enter in 
the cross-section for the inclusive photon Eq. (|15p and semi-inclusive photon-hadron Eq. (|12p. and single inclusive 
hadron production Eq. ([6]) are different. 



III. DISCUSSION AND PREDICTIONS 



We first present our results for the total charged hadron multiplicity distribution in p-l-Pb collisions at ^fS = 5.02 
TeV. The details of the calculation can be found in Ref. Our results is consistent with predictions given in 

Ref. 15] at 4.4 TeV for minimum-bias collisions. Here, we quantify the theoretical errors and provide predictions 
at 5.02 TeV at different centralities. The main input in the kT factorization for calculating the multiplicity is the 
dipole forward amplitude which we take the b-CGC saturation model [3l|. This model explicitly depends on the 
impact-parameter and approximately incorporates all known features of small-x physics and describes the small-x 
data at HERA including diffractive data [HH, as well as RHIC and the LHC data at small-x [Ull^, HI HI) ■ We 
have only two free parameters here, the mini-jet mass and the over-all normalization factor in Eq. ^ which are 
fixed at lower energy^ [l5| . Unfortunately, these two parameters cannot be uniquely fixed, because the error bars in 
experimental data for multiplicity at lower energy is rather large, and moreover, there is a correlation between these 
two parameters. These uncertainties are carefully quantified and shown in Fig. [1] The same setup provides excellent 
description of charged hadron multiplicity distribution in p-l-A collisions at RHIC at different centralities over a wide 
range of rapidity and also in p-l-p and A-|- A collisions at the LHC and lower energies jl5|, [l^l . 

In Fig. [T] (curves labeled by b-CGC), we show the charged hadron multiplicity in p-|-Pb collisions at 5.02 TeV for 
various centralities — 20%, 20 — 40%, 40 — 60%, 60 — 80% and minimum-bias, obtained from the fcr factorization 
Eq. ([1} using the b-CGC saturation model. The impact-parameter dependence of the saturation model is crucial 
here for defining the centrality of the collisions. In Fig. [l] it is seen that the multiplicity distribution has the biggest 
asymmetry for more central collisions while for peripheral collisions like for example at 60 — 80% centrality cut, the 
system becomes more similar to p-l-p collisions, and this fact is reflected in the total charged hadron multiplicity 
distribution. This effect can be more clearly seen in left panel in Fig. [TJ Note that in all predictions shown in Fig. [TJ 
we assumed a flxed mini-jet mass mjet = 5 MeV for all energies/rapidity and centralities. As we already pointed out, 
the mini-jet mass is related to pre-hadronization/hadronization stage and cannot be obtained from saturation physics, 
and it was fixed via a fit to lower energy minimum-bias data. Unfortunately, for very peripheral collisions where the 
system becomes more like p-fp collisions (or symmetric), this assumption is less reliable. More importantly, one 
should also note that we used the fcy factorization which is only proven for asymmetric collisions like p-t-A collisions 
at small-x, therefore for more peripheral collisions the current CGC prescription is less reliable. 

The upcoming experimental measurements of the charged hadron multiplicity at different centralities at the LHC 
can provide vital information about the impact-parameter dependence of the saturation scale. This can be also 
considered as a non-trivial test of underlying dynamics of gluon saturation, as the system at different centralities 
evolves from an asymmetric to more symmetric form and the total charged hadron multiplicity changes significantly, 
roughly speaking, up to a factor of about two compared to minimum-bias collisions, see Fig. [T] 

The recent LHC data on the multiplicity in A-l-A collisions shows that the energy growth of multiplicities in A-l-A 
collisions is different from the corresponding p-l-p collisions [s^l- At first sight, the different power-law behavior of 
multiplicity in A-l-A and p-|-p collisions seems to be in conflict with the very idea of universality of particle production 
in the CGC/saturation picture. However, one should bear in mind that the fcx factorization has been proven only 
at leading log approximation for a scatterings of a dilute partonic system on a dense one such as p-|-A collisions at 



^ Note that the over-all normalization parameter absorbs three unknown factors, the X-factor, a factor due to parton-hadron duality and 
a factor to relate the interaction area to the average impact-parameter of interaction. 
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b-CGC [p+Pb, 5 .02 Te V] b-CGC [p+Pb, 5 .02 Te V] 




-3-2-10123-3-2-10123 
ri ri 



FIG. 1: Pseudorapidity distribution of the charged particles production in p+Pb coUisions at the LHC x/S = 5.02 TeV at 
various centrahties, in right panel, from top to down: — 20%, 20 — 40%, minimum-bias, 40 — 60%, 60 — 80%. Left panel: for a 
better comparison, we show again theoretical curves for more peripheral centrality bins (the theoretical curves in both panels 
are the same). The theoretical uncertainties of about 5% due to fixing the over-all normalization at RHIC are also shown. 



high-energy [l^ [IJ] . On the contrary, it has been already shown that in the case of dense-dense collisions such as A-|- A 
central collisions it is not valid [33] . In Ref. [2l| , we showed that the gluon-decay cascade in pre-hadronization in the 
so-called MLLA (Modified Leading Logarithmic Approximation) regime which has a different kinematics compared 
to the BFKL type gluon emission, gives rise to extra energy-dependence. Although this should be still considered as 
initial-state effect, it is closely related to the open question of how to accommodate the fragmentation processes in 
the kT factorization formalism. In Ref. we extracted gluon-jet decay cascade from e~^e^ annihilation data and we 
showed that the energy-dependence of about 5*" °^^ due to gluon-decay cascade is exactly what explains the different 
power-law energy-dependence of hadron multiplicities in A-f A compared to p+p collisions at the LHC. This effect is 
more important for A-l-A collisions at high energy where the saturation scale is larger and consequently the average 
transverse momentum of the jet becomes larger than 1 GeV. The MLLA gluon decay cascade enhances the hadron 
multiphcity about 10 - 25% [2l|. 

A comparison of various predictions coming from different saturation models for the charged hadron multiplicity 
in minimum-bias p-f A collisions at 4.4 TeV, can be found in Ref. [34]. The upcoming p-f A collisions at the LHC can 
in principle discriminate between different approaches and put more constrain on saturation models. 

Next we provide our prediction for spectra of the single inclusive hadron and prompt photon production in terms 
of the nuclear modification factor^ RpA defined as follows. 



Ncoii (Pprdrj cPprdr] 



where Ncoii is the number of binary proton-nucleus collisions. We take Ncoii = 6.9 in minimum-bias p-|-Pb collisions 
at -s/S = 5 TeV [s^. In order to compare our RpA predictions with experimental data, one needs to rescale our 
theoretical curves by matching the normalization Ncoii to the experimental value. We will use the NLO MSTW 2008 
PDFs ^3^ and the NLO KKP FFs [13] for neutral pion and charged hadrons. For the photon fragmentation function, 
we will use the full leading log parametrization [sO, [11] . We assume the factorization scale Q in the FFs and the 



^ Upon request, we can also provide our theoretical predictions for the spectra of inclusive hadron production in p+A and p+p collisions 
at the LHC. 
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PDFs to be equal. In order to investigate the uncertainties associated with the choice of the hard-scale Q, we consider 
various cases of Q = 2pt,PttPt/'^ and Q = 2p7p,p7p^p7p/2 for inclusive hadron production Eq. ([5]) and prompt photon 
production Eqs. (|13I14|) . respectively. 

Solving the rcBK equation in the presence of the impact-parameter is still an open problem. However, for 
the minimum-bias analysis considered here the impact-parameter dependence may be less important. Then, the 
initial saturation scale Qos should be considered as an impact-parameter averaged value and it is extracted from the 
minimum-bias data. The main input in the rcBK evolution equation is the initial dipole profile Eq. Q with two 
unknown parameters Qos and 7. The current HERA and RHIC data alone are not enough to uniquely fix the values 
of Qop and 7, namely there exists several parameter sets of {Qop, 7) which all provide a good description of the data^ 
[31 ■ However, the recent LHC data on p-l-p collisions for spectra of the inclusive hadron production provided more 
constrains on the values of Qop and 7. It was recently shown that Qq^ « 0.168 GeV^ with 7 » 1.119 provides a good 
description of the small-x data at the LHC, HERA and RHIC with a proton target IHjESI) see also Ref. [l^. We 
will take these values for the rcBK description of the projectile proton. 

Having known the rcBK description of the proton, now we should determine the rcBK description of a nucleus. 
In the CGC picture, at high-energy or small-x, there is no difference between a proton and nucleus except in their 
saturation scale. Obviously in the case of the nucleus, the impact-parameter dependence is more important. Given 
that the perturbative rcBK evolution equation does not incorporate non-perturbative confinement effect [lo| and the 
b-dependent solution of the BK equation is not yet known, it seems like we are doomed to go beyond the current CGC 
framework and add several model assumptions. Here, we stay in the standard CGC framework without invoking any 
extra ingredient to the model. Our only assumption is that like in the case of proton, the initial saturation scale of a 
nucleus Qoa is an impact-parameter averaged value and can be extracted from other reactions with nuclear targets. 
In our approach, the role of the fluctuations on particle production and all other possible non-perturbative effects 
are effectively incorporated into the average value of Qqa extracted from experimental data. Then we assume that 
the small-x evolution (the BK or JIMWLK equations) will not be altered due to non-perturbative soft physics. For 
the minimum-bias collisions, it is generally assumed that the initial saturation scale of a nucleus with atomic mass 
number A, scales linearly with A^/^ [1, SISIj namely we have 

QIj, = q2^, (17) 

where the parameter c can be fixed via a fit to data. In Ref. ji^], it was shown that New Muon Collaboration's 
(NMC) data at small-x can be described with c sa 0.5. This is also consistent with the fact that RHIC inclusive 
hadron production data in minimum-bias deuteron-gold (d-f A) collisions and DIS data for heavy nuclear targets 
can be described with sa 3 4 Ql^ [H, Hi], [H, HJ. However, one should bear in mind that RHIC and DIS 
data on heavy nuclear target are very limited in number of data points with rather large errors. Most importantly, 
RHIC data in d-|-A collisions at forward rapidities and nuclear DIS data at small-x are limited in kinematics to low 
transverse momenta (and low virtualities) , about pr < ^ GeV. At high px, we expect that — ?> 1. Assuming that 
in high-energy collisions at high-px [pt > 4 GeV) we are still in the saturation region (which can be the case at the 
LHC), then by making use of the fcr factorization or hybrid formalism, at high-p^ we approximately have, 

where Sa and Sp are effective interaction area in p+A and p+p collisions, respectively. This leads to the same relation 
as in Eq. P7|) . but with c « 1. Therefore, for heavy nuclei (with A « 208), the initial nuclear saturation scale can 
vary within: 

3^4g2p<g2^(a;o = 0.01)<6-7Q2^. (19) 

We recall that the initial saturation scale of proton extracted from HERA data is rather small Qq.^ — 0.168 GeV^. 

Therefore, the upper limit of the initial saturation scale of lead nucleus is about Qqj^{xo = 0.01) < 1.2 GeV^ which 
is consistent with phenomenological saturation models based on RHIC and HERA data, see Ref. 3] and references 
therein. In Eq. (IT9l) . it seems as if we made an assumption at which value of x we will have i?^^ 1. However, 
the assumed value of xq will not affect the inequality given in Eq. (|19p in saturation region. In other words, because 
the speed of the rcBK evolution in rapidity is the same for both proton and nucleus, the inequality given in Eq. (jl9p 



This is partly because there is indeed correlation between parameters of the model, and in order to constrain the model, more exclusive 
data at larger kinematic regions are required. In this sense, the upcoming LHC data on the p+A collisions will be very useful. 
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should be also valid at a;o = 0.01. In previous studies of the nuclear modification at the LHC 0,0,113, in obtaining 
the rcBK equation solutions, only lower limit of the initial saturation scale for nucleus was considered. Here we will 
extend our analysis by considering entire range of Qoa given in Eq. (fT9|) . 

In Fig. [21 we show the nuclear modification Rp'X for inclusive charged hadron production in minimum-bias p+Pb 

collisions at = 5 TeV at different rapidities r] = 0, 2, 4, 6 obtained from the hybrid factorization Eq. ([5]) supple- 
mented with the rcBK evolution equation. The band labeled CGC-rcBK includes uncertainties associated with the 
variation of initial nuclear saturation scale within the range given in Eq. (1191) in obtaining the rcBK equation solu- 
tions, and also the uncertainties coming from various choices of the factorization scale Q = 2pt-,pt,Pt/'^ in Eq. ([5]). 
It was already shown in [l^ that unfortunately with RHIC data alone, one cannot fix the value of a^" in Eq. ([5]), 
and with a reasonable ii'- factor any value of a™ within 0.3^0.2 < a™ < can be consistent with RHIC inclusive 
hadron production data both in p+p and p-|-A collisions. In Fig. [21 we show the effect of various values for the 
strong-coupling a™ in Eq. ([5]), namely in every panel for a given rapidity we also show the results with a™ = and 
0.3 0.2 < a™ < 0.1. In Fig. [2l it is generally seen that inelastic contribution to the hybrid formalism distorts the 
Rp'X flatness with transverse momentum, namely at \ow-pT for more central collisions it leads to more suppressions 
while for more forward collisions it enhances Rp'X at highp-r. This behavior was numerically first observed in Ref. [l6j. 
It is interesting to notice that for more central collisions, let's say at 77 = 0, 2, where the higher-order corrections are 
important, the inelastic contribution becomes more sensitive to the saturation effect and it leads to more suppression 
at low pt- This is more obvious for a bigger initial nuclear saturation scale at low-pr- This is in accordance with the 
fact that inelastic terms are generally more sensitive to the saturation physics. This expectation was first qualitatively 
pointed out in Ref. j27| . 

In order to better quantify the effect of uncertainties due to different initial nuclear saturation within the constrained 
range given in Eq. (IT5)) . in Fig. [21 we show i?^^ obtained from rcBK solutions with different initial nuclear saturation 
scale Qq^ = NQqp with 3 < < 7, and with factorization scale taken Q = pr- In order to demonstrate uncertainties 
associated with various value for the strong-coupling a™ in Eq. ([5]), we also show i?^^ for fixed values a™ = 0, 0.1, 0.2. 
Unfortunately, i?^^ is very sensitive to the initial saturation scale of nucleus, and variation of Qqa in range given 
in Eq. (|T9l) gives rise to rather large uncertainties for RpX- However, if the experimental measurement of Rp\ at 
one rapidity is performed then one can use this data to extract the initial nuclear saturation scale from Figs. (|2l3p . 
Let's assume that experimental value of i?^^ at one rapidity, let's say at = 2, lies between two lines labeled by 
Ni,N2, then for other rapidities, our predictions between two lines with the same labeling Ni,N2 which has the same 
Qoa, should be only considered as true predictions of the CGC. In this way, by knowing i?^^ at one rapidity, one 
can significantly reduces theoretical uncertainties associated with Qqa at other rapidities. Note that our predictions 
here are consistent with earlier studies in pj, 40, 45j. It is generally seen that the CGC predicts more suppression 
at forward rapidities compared to the collinear factorization results Moreover, the small-x evolution washes 

away the Cronin-type peak at low-pT at all rapidities. Therefore, observation of a strong Cronin-type peak in p-f-Pb 
collisions at the LHC, regardless if it is enhancement or suppression, can be potentially a signal of non-CGC physics, 
and it will be difficult to accommodate this feature in the CGC approach at current accuracy*, see also Refs. |47h50| . 

In Figs. p|4p . we show our predictions for the nuclear modification i?^^ of direct photon and single inclusive prompt 
photon production in minimum-bias p-|-Pb collisions at 5 TeV at various rapidities obtained from Eqs. (|13I14|) with 
solutions of the rcBK with different initial saturation scale for nucleus. Similar to the case of inclusive hadron 
production, i?^^ for direct photon production (and inclusive prompt photon production) is also very sensitive to the 
initial saturation scale of nucleus. The band labeled by CGC-rcBK includes uncertainties due to the variation of the 
initial saturation scale of nucleus and different factorization scale Q = 2p^,p'^,p^/2. We also show our predictions 
for with various Qq^ = NQ^p and = 3 7 constrained in Eq. (|19p . Again, the value of A^ can be extracted 

from Rpj^ or Rp\ measurements at one rapidity using the predictions shown in Figs. (|2|3I4|) for different values of 
N. Then at other rapidities the corresponding curves labeled with N should be taken as our prediction. It is seen 
from Figs. (|3l4p that i?^^ of direct photon and inclusive prompt photon are rather similar, although the direct photon 
production is generally more suppressed than the inclusive prompt photon production. This is what one may expect 
in our picture since direct photon cross-section in Eq. ()13p probes the target structure function at lower transverse 
momentum p^ and consequently lower x than the fragmentation part. 

In Fig. [4l at = we also show i?^^ for the inclusive prompt photon production calculated in 147] by lancu- 
Itakura-Munier (IIM) saturation model jSlj. In the IIM model, the saturation is approached from the BFKL region, 
and therefore the small-x evolution encoded in this model is different from the rcBK equation. Note that the IIM 



This was also pointed out in the talks given by Xin-Nian Wang and Adrian Dumitru in pA@LHC workshop, CERN, Geneva, June 2012. 
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FIG. 2: The nuclear modification factor R'jIa inclusive charged hadrons + h production in minimum-bias p-f Pb 
collisions at x/S = 5 TeV at different rapidities = 0, 2, 4, 6 obtained from the hybrid factorization Eq. (|19|l with the solutions 
of the rcBK with different initial saturation scale for nucleus. The band labeled CGC-rcBK includes uncertainties due to the 
variation of the initial saturation scale of nucleus and different factorization scale Q. At every rapidity, we also show the 
results by taking a™ = (only elastic contribution) and 0.3 -r- 0.2 < a™ < 0.1. The lines labeled by a number TV are the 
results with a fixed hard-scale Q = Pt and a fixed saturation scale Q^a ~ NQop with A'^ = 3 7 constrained in Eq. (|19() and 
Qlp = 0.168 GeV^ 
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FIG. 3: The nuclear modification factor i?^^ for direct photon production in minimum-bias p-|-Pb collisions at \fS = 5 TeV at 
different rapidities r\~f — 0, 2, 4, 6 obtained from Eq. 1)14^ with the solutions of the rcBK with different initial saturation scale 
for nucleus. The band labeled CGC-rcBK includes uncertainties due to the variation of the initial saturation scale of nucleus 
and different factorization scale Q. Similar to Fig. [2] the lines labeled with a number A*' are the results with a fixed hard-scale 
Q — pip and a fixed saturation scale Qqa ~ ^Qop with A'^ = 3 7 constrained in Eq. (|19p and Qop ~ 0.168 GeV^. 



saturation model also provide a good description of HERA data [5i[ . Nevertheless, it is seen from Fig. |3]that depending 
on the value of Qqa, the rcBK and the JIM saturation model provide different suppression for inclusive prompt photon 
production in minimum-bias p-|-A collisions at the LHC^. Note that similar to the inclusive hadron production, the 
CGC approach generally gives larger suppression for inclusive prompt photon i?^^ at forward rapidities compared to 
coUinear factorization results (sj]. 

In Fig. [5l we show the effect of different choices for the hard-scale Q appeared in the factorization formulas 
Eqs. (|5ll4p for both inclusive hadron production (right panel) and direct photon production (left panel) at forward 
rapidity rj — rj-y — 2. The initial nuclear saturation scale is fixed for all lines. Note that in the case of inclusive hadron 
production with a™ = (only elastic terms), the hard-scale Q only appears in the PDFs and FFs, and the effect of 
different value for Q is negligible in the Rp'X- However, in the presence of inelastic contribution (a™ 7^ 0), the choice 
of hard-scale Q becomes important and it leads to a sizable effect for the RpX, see Fig. [S] The sensitivity of the 
nuclear modification factor i?^^ and i?^^ to the hard-scale Q (and the strong-coupling a™) clearly indicates that the 
higher-order corrections should be important. Note that for the case of inclusive hadron production, the full NLO 
corrections to the hybrid formalism has been recently calculated [ssj . but yet to be employed for phenomenological 
purpose. 

Some words of caution are in order here. Our formulation is valid for asymmetric collisions when a projectile can 
be treated in the standard collinear approximation while for the target we systematically incorporated the small-x re- 
summation (at the leading twist approximation) effects. Note, however, our reference for or BZ^ is p-l-p collisions 



^ Note that in Ref. |47|| . the coordinate representation of Eq. l llOI I was used. But as it was shown in Ref. [29II . this is equivalent to the 
CGC formulation given in Eq. JTOj, see also Refs. [52ll53(|. 
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FIG. 4: The nuclear modification factor R^^ for single inclusive prompt photon production in minimum-bias p+Pb collisions 

at y/S = 5 TeV at different rapidities r]^ — 0, 2. The descriptions of the band and the solid black lines are the same as in Fig. |3] 
The dashed red line (labeled CGC-IIM) is calculated from the same master Eq. (|10|) but a different saturation model, namely 
the so-called lancu-Itakura-Munier (IIM) saturation model [13, HJl . 




FIG. 5: The nuclear modification factor in minimum-bias p-|-Pb collisions at the LHC \/S — 5 TeV and rj = t]-., = 2 for 
inclusive hadron production (right panel) and direct photon production (left panel) with different hard-scale (or factorization 
scale) Q in Eqs. (I14I5|I . The initial saturation scale for nucleus is fixed for all lines with N = 7. The bands are the same as in 
Figs.lHEl). 



where at mid-rapidity and low transverse momentum the interacting system becomes symmetric, and therefore the 
CGC formulation (both the fc^ factorization and the hybrid formalism) will be less reliable. Moreover, our parameter 
sets for the rcBK equation was obtained from a fit to HERA data at small-x x < 0.01 and for low virtualities 
Q2 e [0.25, 40] GeV^ Q. Therefore, our predictions at high-pT {pt,Pt > 6 7 GeV) should be taken with a grain 
of salt. 

Finally, in Fig. |6l we show our predictions for the azimuthal photon-hadron correlations in minimum-bias p+p and 
p-|-Pb collisions at 5 TeV for two different kinematics, namely in right panel pi^ > pip and in left panel p^ < pip. The 
azimuthal photon-hadron correlation is defined as [TtI . [18| , 

P(A0) = ^ [A0]/ ^ [A(b = T:/2], (20) 

(Pbt ppdpp ppdpp drj-y dij^ d(f> d'^bt Ppdpp ppdpp drj^ drj^ d(j) 
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FIG. 6: The photon-hadron correlation P(A<j!>) defined in Eq. (|20p in minimum-bias (Mini-Bias) p+p and p+Pb collisions 
at \/S = 5 TeV and rjh = rjj = 3 obtained via the rcBK evolution equation with different initial nuclear saturation scale 
QoA = NQqp and A*' = 3, 5, 7 for two different transverse momentum regions > (right) and < Pt (left). 



where the photon-hadron cross-section is given in Eq. (|lip . The function P(A0) gives the probability of the semi- 
inclusive photon-hadron pair production at a certain kinematics and angle A0, triggering the same production with 
the same kinematics at a fixed reference angle Afjjc = tt/2. It is seen that given the transverse momenta of the 
produced photon and hadron, the photon-hadron correlation can have a double or single peak structure. In Ref. [l^ 
it was shown that this feature is related to saturation physics and it is mainly controlled by the ratio p'^/p]- (as can 
be seen in Fig. (5]). We refere the interested readers to Ref. [Ts'l for the details and discussions. Here we only stress 
that the photon-hadron azimuthal correlations in p-l-A and p+p collisions can be considered as an excellent probe 
of the small-x dynamics [13, [13] ■ In Fig. [6l we show the sensitivity of the correlations to various initial saturation 
scale. Although the uncertainties associated to the initial saturation scale changes the strength of the correlation, 
remarkably the main features of the photon-hadron correlations which are intimately connected to saturation physics, 
seem to be robust, namely the suppression of the away-side correlations in p-l-A compared to p+p collisions, and 
appearance of the double or single peak structure. The predictions for the azimuthal angle correlations between the 
produced prompt photon and hadron calculated via the coincidence probability in p+A collisions at y/S = 5 TeV can 



be found in Ref. |18|. 



To summarize, in this letter within the CGC/saturation framework, we provided various predictions for different 
observables for the upcoming p+Pb collisions at ^fS = 5 TeV at various rapidities. The main caveats of the current 
CGC approach are twofold, lack of systematic control of higher order corrections and lack of enough high-quality 
experimental data at small-x to constrain the initial dipole profile of the rcBK evolution equation. We showed that 
the nuclear modification factor for inclusive hadron and direct photon production at the LHC are sensitive to both 
caveats. Unfortunately, the current small-x data on heavy nuclei cannot uniquely fix the dipole parameters and the 
initial nuclear saturation scale Q()a{x — 0.01). This gives rise to rather sizable theoretical uncertainties for RpA- 
We made detailed predictions for different observables in p+Pb collisions at the LHC with various values of Qqa 
(labeled by N) constrained by existing experimental data at small-x. If experimental data at a given rapidity for 
observables considered here are know, one can extract QoA (or N) by confronting our predictions in Figs. (|2I3I4|6|) with 
experimental data, then at other rapidities our results labeled with the corresponding TV should be only considered 
as true CGC predictions. In this way, despite our rather large theoretical uncertainties, one can still test the main 
dynamics of the CGC/saturation at the LHC. 
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Note added 29 Nov. 2012: Comparing our predictions with ALICE data 

After we submitted our paper to archive, the ALICE collaboration released the first p+Pb data from a short pilot 
run performed in preparation for the full p+Pb physics run scheduled for the beginning of 2013. Here, we compare 
some of our predictions already presented in this paper with the ALICE preliminary data [56l.[57j. In Fig. [71 we show 
our results for the charged-particle pseudorapidity density in non-single diffractive p+Pb collisions at = 5.02 
TeV. In order to compare with ALICE data, we have also accounted for the boost of the 77 = lab frame by adding 
a rapidity shift of Ay = —0.465. The details of calculation can be found in Sees. II, III. We recall that the main 
theoretical uncertainties in our approach is due to our freedom to choose the value of mini-jet mass rrijet appeared 
in Eq. Unfortunately, RHIC data alone is not enough to uniquely fix the value of rrijet ■ In Fig. [71 we show 
our results with different values of mjef All assumed values of mini-jet within mjet = 0.001 0.03 GeV give a 
good description of RHIC data for charged hadron multiplicity. In Fig.[7l (right panel), we show our prediction with 
nijet = 0.03 GeV which was available before the ALICE data and it is in agreement with data within about (or less 
than) 10% errors. However, it appears that rrijet ~ 5 MeV gives the best description of the ALICE data with less 
than 4% errors (see Fig. [7] left panel). This value is remarkably similar to current quark masses. In Fig. [l] we also 
assumed a fixed rujet = 5 MeV for all different centralities. We recall that using the kx factorization, one needs to 

rewrite the rapidity y distribution in terms of pseudorapidity via y{h) — i log ^'^"^^ with the Jacobian 

Y cosh^ ?7-|-^i2 — sinh ?7 

of rapidity-pseudorapidity transformation h = dy/drj. The scale /i is determined by the typical mini-jet mass and 
its transverse momentum. Note that different CGC calculations based on the fc^ factorization but using different 
saturation models (KLN (5^, IP-Sat [11], rcBK [i^j and b-CGC in this paper) have employed different definition for 
the scale ^: 



2 

2 _ ™jet 



= ^ = 0-24- - 0.0035, [iV,.,.,P,-l]^ ^, 
(0.13 + 0.32(V^/lTeV)oii5)2 ^ 



rr? 

^2 ^ _jet^ with rnjei = 0.2 + 0.4 GeV. Ref. [ii]: IP-Sat. 

Pt 

= = . Ref. rcBK. 

(0.13 + 0.32(V^/lTeV)oii5)2 ^ 

2 

2 ^^jet 

jJL = — with mjet = TOcurront quark ~ 0.001 + 0.01 GeV. This paper: b-CGC. 

Pt 

Unfortunately, the value of rrijet is interconnected with both soft and hard physics and its true value cannot be 
determined at the current approximation. We recall also that the value of mini-jet mass changes the over-all prefactor 
behind the Ut factorization indicating that rrijet may mimic some higher order corrections. Nevertheless, our freedom 
to chose different values for /x or rrijet may bring uncertainties as large as 5 + 15% at the LHC. Note also that 
currently all the CGC predictions for charged hadron multiplicity have been based on fixed-coupling /cT-fa-ctorization 



15 




FIG. 7: Pseudorapidity distribution of the charged particles production in minimum-bias (Mini-bias) p-|-A colhsions at the LHC 
VS = 5.02 TeV. The theoretical curves labeled by b-CGC is based on leading log fct-factorization formalism and the b-CGC 
saturation model. In the right and the left panel we show the results with mini-jet mass rrij^t = 30 MeV and rrijet = 1,5, 10 
MeV respectively. The theoretical uncertainties of about 5% due to fixing the over-all normalization at RHIC are not shown. 
The experimental data are from the ALICE (preliminary) 56], PHOBOS and BRAHMS collaborations [s^ . 
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FIG. 8: The nuclear modification factor Rp\ for inclusive charged hadrons production in minimum-bias p+Pb collisions at 
y/S = 5.02 TeV at 77 = obtained with the solutions of the rcBK with average initial saturation scale for nucleus Qqa = 
iV 0.168 GeV^ with TV = 4 6 (grey area) and iV = 5 (black curve) constrained in Eq. (| 19]) . The theory curves are taken from 
Fig. (2] The preliminary (prelim.) experimental data are from the ALICE collaboration [57|]. The dashed lines correspond to 
maximum systematics errors of experimental data points. 



formula while incorporating the running-coupling phenomenologically. However, recently, running coupling corrections 
for the lowest-order gluon production in high energy hadronic and nuclear scatterings was calculated [60] and it 
was conjectured how running coupling corrections may enter the full fixed-coupling fcy-factorization formula. This 
formulation is yet to be explored for future phenomenological applications. 

Finally in Fig. [51 we compare our predictions shown in Fig. [2] with preliminary ALICE data for the nuclear 
modification factor i?^^ of inclusive charged hadrons production in minimum-bias p+Pb collisions at \/S = 5.02 TeV 
at ?7 = 0. The ALICE data are in excellent agreement with our predictions (shown in Fig. [2]) with the solutions of the 
rcBK evolution equation with average initial saturation scale for nucleus = TV 0.168 GeV^ with N ^ 5 constrained 
in Eq. (ITO)) . It is remarkable that the preferred value of N corresponds to the average value of Qqa extracted from 
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other reactions given in Eq. p9|) . Unfortunately, ALICE data points have rather large systematics errors, nevertheless, 
it is seen that ALICE data already put a strong extra constrain in Eq. and prefers value within TV « 4 6 with 
a™ w (shown in Fig. [2] with grey area). Although, a bigger N with larger a™ also cannot be ruled out at the 
moment. As we already stressed the scale of a™ cannot be determined with the current approximation and it requires 
a full NNLO calculation which is not yet available. Therefore, our freedom in choosing the value of a"' in hybrid 
factorization formula Eq. ([T9l) brings rather large uncertainties. One of the most remarkable feature of the ALICE 
data for i?^^ is the fact that data seem to rule out any (or strong) Cronin-type peak. Although, the experimental 
error bars are still rather large to make any firm conclusion, but if this feature of data remains intact with more precise 
experimental data, this can be considered as another important evidence in favor of the CGC approach as this feature 
was already predicted, see discussion in Sec III. The upcoming measurement of Rp\ at forward rapidities at the LHC 
can be considered as a crucial further test of the CGC approach and provide undoubtedly valuable information about 
the saturation dynamics and the true value of the initial saturation scale of heavy nucleus Qqa- 



